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Introduction 

Characterizing dust plumes on the moon's surface during a rocket landing is imperative to the 
success of future operations on the moon or any other celestial body with a dusty or soil surface 
(including cold surfaces covered by frozen gas ice crystals, such as the moons of the outer plan- 
ets). The most practical method of characterizing the dust clouds is to analyze video or still cam- 
era images of the dust illuminated by the sun or on-board light sources (such as lasers). The 
method described below was used to characterize the dust plumes from the Apollo 12 landing. 

In this context, the histogram matching method (HMM) is an image processing technique for de- 
termining dust optical density in Apollo landing videos. The software implementation of HMM 
creates a greyscale image histogram and calculates the histogram mean and standard deviation, 
which is then used to match dusty and clear images for the purpose of estimating an effective 
optical density and optical extinction coefficient a. A dust thickness model, based on the tilt of 
the camera and increasing height of the dust layer towards the top of the image, is used to ac- 
count for the distance light travels through the dust. 

Previous methods relied on comparing specific features in clear vs. dusty images, which severely 
limited ability to analyze video frames. This method compares the statistical nature of a clear im- 
age to the statistical nature of a dusty image, assuming that the average scene's description (as 
characterized by an image histogram) due to surface reflectance and sun angle is invariant 
throughout the frame sequence. This assumption fails when shadows show up on the scene, 
which is evident in the last 20 s of the landing descent. In the last 20 s, the error minimization of 
the histogram matching is by-passed and the matching is done manually by visually comparing 
images. 

Regolith erosion rate and total mass eroded due to a vehicle landing on the lunar surface can be 
computed by a technique similar to that used by the National Weather Service for estimating 
rainfall rate from radar reflectivity [1]. With this proposed method, soil erosion rate is similar to 
rainfall deposition rate, while optical extinction [2] is similar to weather radar reflectivity. Radar 
meteorology relies on a predetermined average rainfall drop size distribution in order to extract 
rainfall rate from radar measurements. Similarly, the surface regolith erosion rate due to plume 
impingement of a landing vehicle [3] depends on knowledge of the soil’s particle size distribu- 
tion [4]. In the case of Apollo landers, optical extinction can be estimated by image analysis of 
the cockpit camera landing videos. Using this approach, the computed plume induced soil ero- 
sion rate is: 
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where pi is the bulk density of lunar soil (3100 [kg m' ]). The second moment M 2 is related to 
the measured optical extinction by, a = {tiJA) Qe M 2 , where Qe is the scattering efficiency factor 
for extinction [5], The velocity term, v{t, D) is determined from computation fluid dynamics 
(CFD) analysis of the Apollo LM descent engine, as well as particle trajectory analysis based on 
the CFD simulations [6], See Appendix A for details on modeling v{t, D). S(D) is the normal- 
ized soil size distribution of Apollo soil samples returned from the lunar surface, as measured by 
a Fine Particle Analyzer, sieves, or something equivalent (see Appendix B). An empirical rela- 
tionship is established between shear stress as determined by CFD simulation and erosion rate 
estimated by optical extinction measurements for the case of the Apollo 12 LM. 


Histogram Matching Method 

The output of the HMM algorithm is a modified image, where “dust” has been added (mode 1) 
or removed (mode 0). The output image pixel p'^ is computed from the input image pixel p.j : 
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where a is the optical extinction coefficient of the dust and po is a fitting parameter associated 
with the dynamic range of the image (ideally po = 255 for an 8-bit image). The factor of two in 
the extinction term is the result of light reflecting off of the lunar surface back to the camera. 

The exponent term without the factor of two corresponds to light scattered back to the camera 
from the dust cloud. The distance x in Equation (2) is the effective optical-dust path length model 
along the camera view ray through the dust cloud, corresponding to each ij image pixel in the 
image. It is equal to the physical path length xd of the dust for r < uq, where r is the radial dis- 
tance from the engine nozzle centerline and qq is a parameter. For r > uq, the effective path 
length is xd scaled by the radial dispersion factor: 
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where F is the focal length of the camera (F = 10 mm), d is the pixel width (J « 15 /an), 0 is the 
dust angle relative to horizontal (^« 3°), the camera angle relative to vertical {<f>= 33°), Ic = 
1.2 m is the camera offset distance from the nozzle center line, and h is height of the LM above 
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the surface. The variable q is the vertical distance in the image in pixel units from the ij pixel to 
a horizontal centerline in the rotated camera view: 

q = (j- j N)cos - (i- }M)sm<^ , (5) 

where is the camera rotation angle about the camera axis ((^ « -33°), i is the horizontal pixel 
index, J is the vertical pixel index, N is the total number of horizontal pixels, and M is the total 
number of vertical pixels. 

The left side of Figure 1 shows frame F3077 {h = 34 m) of the cockpit video camera. The right 
side of Figure 1 displays a map of the same field of view for this frame, showing the effective 
optical-dust path length model. Equation (3), which is based on the tilt of the camera and increas- 
ing depth of the dust layer towards the top of the image. The horizontal and vertical axes of the 
plot are in pixel units. The contours are graded in increments of 0.3 m, starting with the mini- 
mum, X = 0.3 m at the bottom (purple in the on-line version), to a maximum, x = 3.0 m at the top 
(red in the on-line version). 

Figure 2 is a similar image set, occurring 46.6 s later at an LM altitude of =1 1 m. The contours 
are graded in increments of 0.017 m, starting with the minimum, x = 0.1 m at the bottom (purple 
in the on-line version), to a maximum, x = 0.25 m near the center (red in the on-line version). 
Note that the video frame numbers F3077 (Figure 1) and F3543 (Figure 2) correspond to a con- 
stant frame rate of 10 fps. 

The HMM algorithm processes two input images, pixel by pixel. The first input image is repre- 
sented by pixel py as shown by Equation (2). A reference image is represented by qij. For mode 
= 0, The HHM algorithm applies the transformation described by Equation (2) to the input image 
Pij (frame with dust), creating an output image p 'ij (artificially removed dust), as shown in Figure 
3. The reference image qy (no dust) is then compared to p 'ij and by matching the average and 
standard deviation of the their histograms, the parameters a, po, and ao are found. Figure 4 
shows a similar example for mode = 1 . 
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Fig. 1. Dust depth model: (left) video camera frame (F3077), with LM altitude A = 34 m; (right) 
effective camera dust length x with radius uq = 46 m. 



Fig. 2. Dust depth model: (left) video camera frame (F3543), with LM altitude A = 1 1 m; (right) 

effective dust depth with radius ao = 6.5 m. 
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Fig. 3. HMM algorithm example for mode = 0: (lower right) image py; (lower left) output image p y; (upper left) reference image qy ; 

(upper right) histograms for the three images. 




qij ; (upper right) histograms for the three images. 




The extinction coefficient a can be estimated from the Apollo videos by adding dust to a clear 
image (before dust appears) and comparing to the dusty image of interest (see Figure 5). Reiter- 
ating from the last section, a = kIA Qe M 2 , where Qe is the scattering efficiency, assumed to be 
equal to 2. The second moment in Equation (1) is the product of the measured extinction coeffi- 
cient and 2! K The total mass rejected (total mass displaced) is Equation (1) integrated over ve- 
hicle descent time and over the area where soil is eroded: 
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where «o(0 is the radius on the surface, referenced to the engine nozzle centerline, where erosion 
is taking place. The assumption inherent in Equation (6) is that erosion is uniform over a circle 
of radius ao{f) and zero outside of that circle. An estimate of «o(0 for the Apollo 12 EM is 
shown in Figure 5, which is an output of the HMM algorithm. Note that in this and all previous 
discussions of S{D) and its moments, it has been assumed that the particle size distribution is 
homogeneous over the extent of measurement, i.e., within a circle of radius at){f), and all tem- 
poral effects due to engine LM altitude and thrust occur instantaneously over this spatial extent. 
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Fig. 5. Apollo 12 optical extinction estimate using histogram matching method. The time inter- 
vals correspond to the voice callouts of LM pilot, Alan Bean. 


The erosion rate should actually vary with the state of the gas flowing across the soil, including 
its shear stress, rarefaction, and turbulence; it should also vary with saltation, including the 
downward flux of larger particles that are too heavy to be carried away by the gas as well as 


-7- 


smaller particles that are scattered back down from the entrained cloud via particle collisions. 
Exam in ation of the sandblasting effects on Surveyor III have shown that the downward flux of 
scattered particles is significant [7] and discrete element computer simulations show the im- 
portant but largely unexplored role of mid-flight particle scattering in enhancing erosion rate [8], 
Influence of the gas upon erosion rate should be greatest in an annular region around the vehicle 
[9] while the influence of saltation may be greater in another annulus with larger radius since 
particles travel downrange before striking the surface. Thus the net erosion rate may be some- 
what more uniform and spread over a broader region than if gas effects alone are consid- 
ered. The details of erosion physics — especially in lunar rocket exhaust conditions — are not yet 
well understood, so a constant erosion over a fi nite area assumed, as it is the simplest model and 
therefore a sensible first step. 

Conclusions 

A method for estimating lunar soil erosion rate due to plume impingement of the Apollo 12 Lu- 
nar Module Intrepid during its descent to the lunar surface has been presented. The observables 
are optical extinction and particle size distributions of soil samples returned from the lunar sur- 
face. The optical extinction is measured between the camera mounted inside of the cockpit win- 
dow and the lunar surface during landing. CFD analysis of the Apollo LM descent engine, as 
well as particle trajectory analysis based on the CFD simulations, provides the remainder of the 
necessary data. 



Fig. 6. Soil mass erosion rate m versus Optical extinction a , showing power law fit. 


The CFD simulations provide two key pieces of information: the velocity profile of particles as a 
function of starting distance from the engine nozzle, size of the particle, and height of the lander 
from the surface. Taking an approach similar to the problem of estimating rainfall rate from 
weather radar, Equation (1) was presented as the solution to the problem of estimating soil ero- 


sion rate from optical extinction measurements (see Figure 6). In both cases, the particle veloci- 
ties must be known, as well as the particle size distributions. Quantifying the particle velocity 
function is in some sense the most difficult part of the problem for both lunar plume observations 
and weather radar estimation of rainfall. For this reason, the velocity function is likely the 
source of greatest error. 

The optical extinction model described above assumes from the start that the erosion rate is uni- 
form over an area defined by radius ao, which is a function of lander height. The primary influ- 
ence of optical extinction as measured by a reduction of brightness of the surface and increase in 
brightness of the dust cloud, is the spatial dust geometry. In this model, erosion is uniform over 
but diverges radially from all points on the surface within r < ao. No erosion occurs for r > 
ao. This approach is loosely analogous to methods traditionally employed for estimating rainfall 
rate from radar, where spatial and temporal averaging strategies are necessary to correlate rain- 
fall accumulation rates and radar reflectivity. 

The result of the Apollo 12 shear stress to erosion rate relation shows an approximate 5/2 power 
dependence of erosion rate on shear stress. Theoretical considerations predict a linear relation- 
ship [10]. The difference is either a consequence of effects not considered by the theory, or er- 
rors in the measurements. 
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APPENDIX A: Particle Velocity Function 


Particle trajectory modeling, based on CFD simulations of the Apollo LM engine and the lunar 
environment, yield a particle velocity function which can be described by the empirical fit shown 
in Equation (A-1): 


v{h,D) = 


0 log Z)+Cj log h>CQ 

1 0 (^ 1 + ^2 log/?)tanh(Zi3+Zi4 log/?) 




otherwise 


(A-1) 


where bo = 0.2964, b\ = -0.225, 62 = 0.1954, 63 = 5, 64 = 4.343, cq = 2.212, and c\ = 3.53. Equa- 
tion (A-1) is based on the particle trajectory results shown in Figues A-1 through A-3. For every 
point in {h, D} space a distribution of particle velocities is computed from the particle trajectory 
code using an equivalent Monte Carlo distribution of initial particle trajectory starting points, 
height above the surface and horizontal distance from the engine nozzle centerline. 

Because of the distribution of particle velocities for a given h and D, the maximum values origi- 
nate near the engine nozzle. Since the erosion area within r<ao (where r is the radial distance 
from the engine nozzle centerline) is much greater for slower velocities, then it reasonable to ex- 
pect rj to be much smaller than 1 . 
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Fig. A-1. A = 45 m case: 
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Case Cl (/i = 6.1 m): xq — 0.03 m, D = \ //m 
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Fig. A-2. h = 6m case. 
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NoCrater (h = 0.76 m): rro = 0.03 m, D = 1 //m 


NoCrater (h = 0.76 m): xq = 0.03 m, D = 10 //m 



Horizontal Range, y [m] Horizontal Range, y [m] 

NoCrater (h = 0.76 m): xq = 0.03 m, D = 100 /mi NoCrater (h = 0.76 m): xq = 0.03 m, D = 1000 /mi 




Fig. A-3. h= l.5m case. 


APPENDIX B: Particle Size Distribution 

The particle size fraction of the lunar soil at the Apollo 12 site is modeled as a combination of 
two power law functions by fitting Apollo 1 1 and 17 soil sample data, as well as JSC- la simulant 
[ 10 ]: 




(B-1) 


where S^(D) = (D/D^y‘‘' , S 2 (D) = (D/D^y^^,and w(D) = {d/ +\. The fitting constants in Equa- 
tion (B-19) are Di = 4.090x10’’ [m], Bi = 1.8, D 2 = 9.507x10’^ [m], B 2 = 5.6, andD^ = 2.5x10'^ 
[m]. The denominator of Equation (B-1) is equal to 4.784x10'^. Note that all units are kept in 
meters even though the numbers are more aesthetically pleasing in micrometers. The reason for 
doing this is to minimize confusion in the integrals involving S(D). 
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